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A HYBRID HIGH-ORDER METHOD FOR THE CAHN-HILLIARD 
PROBLEM IN MIXED FORM* 
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Abstract. In this work we develop a fully implicit Hybrid High-Order algorithm for the Cahn- 
Hilliard problem in mixed form. The space discretization hinges on local reconstruction operators 
from hybrid polynomial unknowns at elements and faces. The proposed method has several ad¬ 
vantageous features: (i) It supports fairly general meshes possibly containing polyhedral elements 
and nonmatching interfaces; (ii) it allows arbitrary approximation orders; (iii) it has a moderate 
computational cost thanks to the possibility of locally eliminating element-based unknowns by static 
condensation. We perform a detailed stability and convergence study, proving optimal convergence 
rates in energy-like norms. Numerical validation is also provided using some of the most common 
tests in the literature. 
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1. Introduction. Let c d e {2,3}, denote a bounded connected convex 
polyhedral domain with boundary dO and outward normal n, and let tp > 0. The 
Cahn-Hilliard problem, originally introduced in [11, 10] to model phase separation in 
a binary alloy, consists in finding the order-parameter c : H x [0, tp] R and chemical 
potential w : 0. x [0, tp] ^ M such that 


(la) 

dtc — Aw = 0 

in H X (0,fp], 

(lb) 

w = $Xc) — 7 ^Ac 

in 0 X (0,fp], 

(Ic) 

II 

o 

in H, 

(Id) 

dnC = dnW = 0 

on dCl X (0,ip], 


where cq e n Lq{0) such that dnCo = 0 on dO denotes the initial datum, 7 > 0 

the interface parameter (usually taking small values), and $ the free-energy such that 

(2) $(c):=l(l-cY. 

Relevant extensions of problem (1) (not considered here) include the introduction of 
a flow which reqnires, in particular, to add a convective term in (la); cf., e.g., [29, 5, 
7, 8, 31, 30]. 

The discretization of the Cahn-Hilliard equation (1) has been considered in several 
works. Different aspects of standard finite element schemes have been studied, e.g., 
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in [22, 21, 14]; cf. also the references therein. Mixed finite elements are considered 
in [24] . In [35] , the authors study a nonconforming method based on shape func¬ 
tions for the fourth-order primal problem obtained by plugging (lb) into (la). Dis¬ 
continuous Galerkin (dG) methods have also received extensive attention. We can cite 
here [36] , where a local dG method is proposed for a Gahn-Hilliard system modelling 
multi-component alloys, and a stability analysis is carried out; [23], where optimal 
error estimates are proved for a dG discretization of the Gahn-Hilliard problem in 
primal form; [30], which contains optimal error estimates for a dG method based on 
the mixed formulation of the problem including a convection term; [26], where a 
multi-grid approach is proposed for the solution of the systems of algebraic equations 
arising from a dG discretization of the Gahn-Hilliard equation. In all of the above 
references, standard meshes are considered. General polygonal meshes in dimension 
d = 2, on the other hand, are supported by the recently proposed C^-conforming Vir¬ 
tual Element (VE) method of [4] for the problem in primal formulation; cf. also [6] for 
VE methods with arbitrary regularity. Therein, the convergence analysis is carried 
out under the assumption that the discrete order-parameter satisfies a C'°(L“)-like a 
priori bound. 

In this work, we develop and analyze a fully implicit Hybrid High-Order (HHO) 
algorithm for problem (1) where the space discretization is based on the HHO(fc -I-1) 
variation proposed in [12] of the method of [19]. The method hinges on hybrid degrees 
of freedom (DOFs) located at mesh elements and faces that are polynomials of degree 
{k + 1) and fc, respectively. The nonlinear term in (lb) is discretized by means of 
element unknowns only. For the second-order diffusive operators in (la) and (lb), on 
the other hand, we rely on two key ingredients devised locally inside each element: 
(i) A potential reconstruction obtained from the solution of (small) Neumann problems 
and (ii) a stabilization term penalizing the lowest-order part of the difference between 
element- and face-based unknowns. See also [13, 34, 33] for related methods for second- 
order linear diffusion operators, each displaying a set of distinctive features. The 
global discrete problem is then obtained by a standard element-by-element assembly 
procedure. When using a first-order (Newton-like) algorithm to solve the resulting 
system of nonlinear algebraic equations, element-based unknowns can be statically 
condensed. As a result, the only globally coupled unknowns in the linear subproblems 
are discontinuous polynomials of degree k on the mesh skeleton for both the order- 
parameter and the chemical potential. With a backward Euler scheme to march in 
time, the C°(iJ^)-like error on the order-parameter and the L^{H^)-like error on the 
chemical potential are proved to optimally converge as (/i^+^ + t) (with h and r 
denoting, respectively, the spatial and temporal mesh sizes) provided the solution has 
sufficient regularity. 

The proposed method has several advantageous features: (i) It supports general 
meshes possibly including polyhedral elements and nonmatching interfaces (resulting, 
e.g., from nonconforming mesh refinement); (ii) it allows one to increase the spatial 
approximation order to accelerate convergence in the presence of (locally) regular so¬ 
lutions; (iii) it is (relatively) inexpensive. When d = 2, e.g., the number of globally 
coupled spatial unknowns for our method scales as 2 card(J^?i)(fc 4- 1) (with card(J^;i) 
denoting the number of mesh faces) as opposed to card(7)i)(fc-(-3)(fcH-2) (with card(7fi) 
denoting the number of mesh elements) for a mixed dG method delivering the same 
order of convergence (i.e., based on broken polynomials of degree fc-(-1). Additionally, 
thanks to the underlying fully discontinuous polynomial spaces, the proposed method 
can accomodate abrupt variations of the unknowns in the vicinity of the interface 
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between phases. 

Our analysis adapts the techniques originally developed in [30] in the context of dG 
methods. Therein, the treatment of the nonlinear term in (lb) hinges on C'^-in-time 
a priori estimates for various norms and seminorms of the discrete order-parameter. 
Instrumental in proving these estimates are discrete functional analysis results, includ¬ 
ing discrete versions of Agmon’s and Gagliardo-Nirenberg-Poincare’s inequalities for 
broken polynomial functions on quasi-uniform matching simplicial meshes. Adapting 
these tools to hybrid polynomial spaces on general meshes entails several new ideas. 
A first key point consists in dehning appropriate discrete counterparts of the Laplace 
and Green’s operators. To this end, we rely on a suitably tailored L^-like hybrid inner 
product which guarantees stability estimates for the former and optimal approxima¬ 
tion properties for the latter. Another key point consists in replacing the standard 
nodal interpolator used in the proofs of [30, Lemmas 2.2 and 2.3] by the L^-orthogonal 
projector which, unlike the former, is naturally defined for meshes containing polyhe¬ 
dral elements. We show that this replacement is possible thanks to the W^’^-stability 
and approximation properties of the L^-orthogonal projector on broken polynomial 
spaces recently presented in a unified setting in [15]; cf. also the references therein 
for previous results on this subject. 

The material is organized as follows: In Section 2 we introduce the notation for 
space and time meshes and recall some key results on broken polynomial spaces; in 
Section 3 we introduce hybrid polynomial spaces and local reconstructions, and state 
the discrete problem; in Section 4 we carry out the stability analysis of the method, 
while the convergence analysis is detailed in Section 5; Section 6 contains an extensive 
numerical validation of the proposed algorithm; finally, in Appendix A we give proofs 
of the discrete functional analysis results used to derive stability bounds and error 
estimates. 

2. Discrete setting. In this section we introduce the discrete setting and recall 
some basic results on broken polynomial spaces. 

2.1. Space and time meshes. We recall here the notion of admissible spatial 
mesh sequence from [17, Ghapter 1]. For the sake of simplicity, we will systematically 
use the term polyhedral also when d = 2. Denote by "H c Kj a countable set of 
spatial meshsizes having 0 as its unique accumulation point. We consider h-refined 
mesh sequences {Th)he'H where, for all h e H, Th is a, hnite collection of nonempty 
disjoint open polyhedral elements T of boundary dT such that D = UreTh ^ 
h = max^eTfe with hx standing for the diameter of the element T. 

A face F is defined as a planar closed connected subset of D with positive (d—1)- 
dimensional Hausdorff measure and such that(i) either there exist Tj,r 2 s Th such 
that F c dTi n dT 2 and F is called an interface or (ii) there exists T e Th such 
that F d dT n dVl and F is called a boundary face. Mesh faces are collected in 
the set Th, and the diameter of a face F e Th is denoted hy hp. For all T e 7/t, 
Tt '■= {A ^ Th I F d dT} denotes the set of faces lying on dT and, for all F e Tp, 
riTF is the unit normal to F pointing out of T. Symmetrically, for all F e Th, we 
denote by Tf the set of one (if F e Tj^) or two (if F e T}^^) elements sharing F. 
Assumption I (Admissible spatial mesh sequence). We assume that, for all h e T-L, 
Th admits a matching simplicial submesh Th and there exists a real number g > 0 
independent of h such that, for all h £ %, the following properties hold:(i) Shape 
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regularity: For all simplex S e %h of diameter hg and inradius rg, ghg ^ rg; 
(ii) contact-regularity: For allT e 7/t, and all S such that S ^ T, gh^ ^ hg. 

To discretize in time, we consider a uniform partition (t"')osSnssA of the time interval 
[Ojtp] with = 0, = t-p and t" — = r for all 1 ^ n ^ TV (the analysis can 

be adapted to nonuniform partitions). For any sufficiently regular function of time ip 
taking values in a vector space V, we denote by v?" G F its value at discrete time t", 
and we introduce the backward differencing operator 6t such that, for all 1 ^ n ^ 


( 3 ) 




F - F 


n—1 


G V. 


In what follows, we often abbreviate by a < 6 the inequality a ^ Cb with a and b 
positive real numbers and C > 0 generic constant independent of both the meshsize 
h and the time step r (named constants are used in the statements for the sake of 
easy consultation). Also, for a subset X c 11, we denote by (•, ■)x and H-Hx the usual 
L^(A)-inner product and norm, with the convention that we omit the index if A = fl. 
The same notation is used for the vector-valued space L^(A)‘^. 

2.2. Basic results on broken polynomial spaces. The proposed method 
is based on local polynomial spaces on mesh elements and faces. Let an integer 
I ^ 0 he fixed. Let U he a subset of Hu the affine space spanned by U, du its 
dimension, and assume that U has a non-empty interior in Hu- We denote by P^(C/) 
the space spanned by dcz-variate polynomials on Hu of total degree I, and by 
the L^-orthogonal projector onto this space. In the following sections, the set U will 
represent a mesh element or face. The space of broken polynomial functions on Th of 
degree I is denoted by P*(7/t), and is the corresponding L^-orthogonal projector. 

We next recall some functional analysis results on polynomial spaces. The following 
discrete trace and inverse inequalities are proved in [17, Chapter 1] (cf. in particular 
Lemmas 1.44 and 1.46): There is C > 0 independent of h such that, for all T G Th, 
and all Vu G P*(T), 


(4) 

u (f ^ u T yFeTr, 

and 


(5) 

IIVuIIt < Chf^Mr- 


We will also need the following local direct and reverse Lebesgue embeddings (cf. [15, 
Lemma 5.1]): There is C > 0 independent of h such that, for all T G Th, all q,p G 
[1, -l-oo], 

(6) VugP^(T), C ^ ^ (7||u||i5(r). 

The proof of the following results for the local L^-orthogonal projector can be found 
in [15, Appendix A.2]. For an open set U of s G N and p G [1, -l-oo], we define the 
seminorm |•|rv=.p(^7) as follows: For all v G 

\v\wo-p(U) '■= Xj II^“^IIlp((7)i 

, \ol\^i =s 
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where \a\ii '.= ax + ■■■ + ad and . For s = 0, we recover the usual 

Lebesgue spaces LP{U). The L^-orthogonal projector is VF®’^’-stable and has optimal 
bF®’^-approximation properties: There is C > 0 independent of h such that, for all 
T e Th, all s e {0,..., Z -I-1}, all p e [1, -l-ao], and all v e it holds, 

(7) 

and, for all m G {0,..., s}, 

( 8 ) \v — 7r^u|tym,p(T) + h^\v — 7ryn|vym,p(jr,j,) ^ Chfp ^\v\w'’,p{t), 

where denotes the set of functions that belong to bF'"’P(F) for all F e Ft- 

Finally, there is C > 0 independent of h such that it holds, for all F e Fh, 

(9) \/veH^{F), ||n - < C'/i|u|/i-i(iP). 


In the proofs of Lemmas 6 and 13 below, we will make use of the following global 
inverse inequalities, which require mesh quasi-uniformity. 

Proposition 1 (Global inverse inequalities for Lebesgue norms of broken polynomi¬ 
als). In addition to Assumption 1, we assume that the mesh is quasi-uniform, i.e., 


( 10 ) yreTh, gh^hT- 

Then, for all polynomial degree I ^ 0 and all 1 ^ p ^ ^ -l-oo, it holds 

(11) Vrcft, G P'(7)t), WwhWLiin) ^ Ch'^~p\\wh\\Lp{n), 

with real number C > 0 independent of h. 

Proof. Let Wh e P*(7)t). We start by proving that, for all p G [1, -foo], 

(12) Vw/i G P'(7h), WwhWL-^^n) ^ h~p\\wh\\Lp{n), 

which corresponds to (11) with q = -l-oo. By the local reverse Lebesgue embed¬ 
dings (6), there is C > 0 independent of h such that 

yT e Fh, \\wh\\L^(T) ^ \\'^h\\LP{T) ^ Cp ^\\uJh\\LP(n), 

where we have used the mesh quasi-uniformity assumption (10) to conclude. Inequal¬ 
ity (12) follows observing that ||u'/i||L<»(n) = maxTeTh I|it;i||f=o(t)- Let us now turn to 
the case 1 ^ q < -l-oo. We have 

where the conclusion follows using (12). □ 

3. The Hybrid High-Order method. In this section we define hybrid spaces 
and state the discrete problem. 
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3.1. Hybrid spaces. The discretization of the diffusion operator hinges on the 
HHO method of [12] using polynomials of degree (fc + 1) inside elements and k on 
mesh faces (cf. Remark 15 for further insight on this choice). The global discrete 
space is dehned as 

(13) C/^= ( X X ( X ) • 

VTeTh / / 

The restriction of to an element T e 7)i is denoted by . For a generic collection 
of DOFs in |7j), we use the underlined notation Vh — i{'UT)TeThJ {'’^F)FeJ^h) for 
all T e 77i, we denote by Vrp = {vt, (i’f)f€J^t) restriction to U^. Also, to keep the 
notation compact, we denote by Vh (no underline) the function in P^+^(7)i) such that 

Vh\T = I'T VT e Th¬ 
in what follows, we will also need the zero-average subspace 

Ulo ■■= {vf,eUU{vh,l) = 0}. 

The interpolator : 7J^(D) —> |7^ is such that, for all v e 77^(0), 

(14) LhV ■■= ((7r^+X)r€rh) (TFi^)F€Ffe)- 

We define on U_h the seminorm such that 

(15) hh\\l,h ■= W'^hVhf + \vh\l,h, khli.h ■= 

where denotes the usual broken gradient on H^{Th) and the stabilization bilinear 
form si,h on U\ x is such that 

(16) Si,,,(v^,Zft) := Y, Zl >T-Yi^FivF-VT),F'^{zF- Zt))f- 

TeTh FeJ^T 

Using the stability and approximation properties of the L^-orthogonal projector ex¬ 
pressed by (7)-(8), one can prove that is i7^-stable: 

(17) yveH^in), ||i;||fi(0). 

The following Friedrichs’ inequalities can be proved using the arguments of [15, 
Lemma 7.2], where element DOFs of degree k are considered (cf. also [9, 16] for 
related results using dG norms): For all r e [1, -too) if d = 2, all r e [1, 6] if d = 3, 

(18) ^ lb/tl|i,?i- 

The case r = 2 corresponds to Poincare’s inequality. Finally, to close this section, we 
prove that 1| defines a norm on UX.o- 

Proposition 2 (Norm l|j|i,h). The map l|j|i,/i defines a norm onU^Q. 

Proof. We only have to show that = 0 Vi^ = 0. By (18), ||Th||i,?i 

Vh = 0, i.e., wr = 0 for all T e Th- Plugging this result into the definition (15) of || j|i,/t, 
we get SreTfe SfeFt ^f^II^fIIf = 0; which implies in turn vp = 0 for all F e Th- □ 


A HHO METHOD FOR THE CAHN-HILLIARD PROBLEM 


7 


3.2. Diffusive bilinear form and discrete problem. For all T e Th, we 

define the potential reconstruction operator —> P^+^(T) such that, for all 

Pt s is the unique solution of the following Neumann problem: 

(19) Vz)t = -{VT, Az)t + 2 {VF. Vz-nTF)F Vz e 

FE:J-T 

with closure condition 1 )t = (TT,l)r- It can be proved that, for all v e 

H^{T), denoting by the restriction of the reduction map defined by (14) to 

H\T)^Ut, 

(20) {V{p^F^I^v-v),Vz)t = Q VzeP'=+i(T), 

which expresses the fact that ° Lt) is the elliptic projector onto P''+^(r) (and, 

as such, has optimal approximation properties in xhe diffusive bilinear 

form Gh on x is obtained by element-wise assembly setting 

(21) ah{vf„z^):= ^ + Si,ft(p^, z^). 

Ten 


with stabilization bilinear form si^h defined by (16). Denoting by ll-lla,?! the seminorm 
dehned by ah on {7^, a straightforward adaptation of the arguments used in [19, 
Lemma 4] shows that 

(22) yVf, e ul, \\vji,h ^ \\vh\\a,h ^ hhWm 


which expresses the coercivity and boundedness of ah- Additionally, following the 
arguments in [19, Theorem 8], one can easily prove that the bilinear form ah enjoys 
the following consistency property: For all v e (D) n Lq(D), I 5= 1, such that 

dnV = 0 on dfl, 


(23) 


sup 


a,(J^,zJ + (Ar;,z;,)| < 


Remark 3 (Consistency of o/j). For sufficiently regular solutions (i.e., when I = k + 2), 
equation (23) shows that the consistency error scales as h^'^^. This is a consequence 
of the fact that both the potential reconstruction p^^ (cf. (19)^ and the stabilization 
bilinear form Si,h (cf (16)^ are consistent for exact solutions that are polynomials of 
degree (k -f 1) inside each element. In particular, a key point in Si^h is to penalize 
t%{vf — vt) instead of {vf — vt)- A similar stabilization bilinear form had been in¬ 
dependently suggested in the context of Hybridizable Discontinuous Galerkin methods 
in [32, Remark 1.2.4]- 

The discrete problem reads: For all 1 ^ n ^ IV, find (c)[, w)[) e U_h.o ^ ILh such that 

(24a) ((5*c)[, ph) + ahiwl ,^J = 0 e Ul 

(24b) «,V';,) = ($'(c)[),V'/.) + 7V(c[[,^J 

and c° e U_h,o solves 

(25) ah{cl,(p^) =-{Aco,ph) e 

We note, in passing, that the face DOFs in c° are not needed to initialize the algorithm. 



F. CHAVE, D. A. DI PIETRO, F. MARCHE, F. PIGEONNEAU 


Remark 4 (Static condensation). Problem (24) is a system of nonlinear algebraic 
equations, which can be solved using an iterative algorithm. When first order (Newton¬ 
like) algorithms are used, element-based DOFs can be locally eliminated at each iter¬ 
ation by a standard static condensation procedure. 

4. Stability analysis. In this section we establish some uniform a priori bounds 
on the discrete solution. To this end, we need a discrete counterpart of Agmon’s 
inequality; cf. [3, Lemma 13.2] and also [1, Theorem 3]. We define on U'f the following 
L^-like inner product: 

{Vh.Zyfjo.h ■■= {Vh,Zh) + So,h(u^,Z^), 

se),h{yh,Ztf) ■.= ^ Y, hpi'K'hivF- zt))f, 

TeT/i F^(Ft 

and denote by ||•||o,;l and j-jo.h the norm and seminorm corresponding to the bilinear 
forms and SQ^h, respectively. For further insight on the role of Sg^h, cf- Re¬ 

mark 18. We introduce the discrete Laplace operator iJf : iff Uff such that, for all 
V.h s U_ht 

(27) —{LLhV-h>^h)o,h = ah{vf^,Zf() yzf, e uff, 

and we denote by L’fVf, (no underline) the broken polynomial function in 
obtained from element DOFs in 

Remark 5 (Restriction of iJf to C/^ o ^ UXp)- Whenever e FhUh ^ U-h.o- 
To prove it, it suffices to take Zf^ = ffxii (27) (with xn characteristic function of 
D), and observe that the left-hand side satisfies {ijfvh, zfo.h = while, by 

definition (21) of the bilinear form au, the right-hand side vanishes. In what follows, 
we keep the same notation for the (bijective) restriction of iJf to Uf q —> Uf q. 

The following result, valid for d e {2,3}, will be proved in Appendix A. 

Lemma 6 (Discrete Agmon’s inequality). Assume mesh quasi-uniformity (10) . Then, 
it holds with real number C > 0 independent of h, 

(28) 

We also recall the following discrete Gronwall’s inequality (cf. [28, Lemma 5.1]). 
Lemma 7 (Discrete Gronwall’s inequality). Let two reals 6,G > 0 be given, and, for 
integers n ^ 1, let , 6", and y" denote nonnegative real numbers such that 

N N 

^ 6" ^ 2 G VA^eN*. 

n=l n=l 

Then, if x^^ < 1 for all n, letting c" := (1 — x"^)~^, it holds 

N / N \ 

(29) a^ + d^fe’^^exp xG VA^ e N^ 

n=l V n=l / 


We are now ready to prove the a priori bounds. 
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Lemma 8 (Uniform a priori bounds). Under the assumptions of Lemma 6, and 
further assuming that t ^ L for a given real number L > 0 independent of h and of r 
(but depending on 7 ^ j and sufficiently small, there is a real number C > 0 independent 
of h and t such that 


max 

IsSnsSA 



N 

+ (‘J’(cIJ), 1) + + ||cft||L>x.(a) + 2 


^ C. 


Proof. The proof is split into several steps, 
(i) We start by proving that 


(30) 


N 


max 

IsSnsSAT 


IIL/t 


2 

I a.h 


+ mc-,),l))+Y,rM\\l, 

n= 1 


< 1 . 


Subtracting (24b) with '0^ = c)) — ^ from (24a) with = ruff, and using the 

fact that, for all r, s e K, $'(r)(r — s) ^ $(r) — $(s) — j(r — s)^, it is inferred, for all 
1 ^ n ^ that 


(31) l^au{cl,cl - cr^) + rIKIIl, + {<^{cl), 1) ^ -\K - cl 


n-l||2 


+ wcr'),i). 


Notice that ($(c)J;),l) > 0 for all 0 ^ n ^ N by definition (2) of $. Making = 
t{cI — cl~^) in (24a) and using the Cauchy-Schwarz and Young’s inequalities, we 
infer that 


(32) 


\\Ch-Ch 


I|2 


^ o\\lLh\\a,h + '^W&h ~ &h \\a,h- 


Additionally, recalling the following formula for the backward Euler scheme: 


(33) 

it holds 


2x{x — y) = x"^ + {x — yY 


(34) 




cr') 


2 (llP/illa,/i 


„n\\2 , II rj 

II /^ L I 


-h Wa.hJ 


Plugging (32) and (34) into (31), we obtain 


l^klWln + ( 7 ^ - 1) kl - crklh + ylKIll. + 2{Hcl), 1) 

< Teller'll.+ 2(<i>(cr^),i). 

Provided r < 27 ^, the bound (30) follows summing the above inequality over 1 ^ n ^ 
N, and using the fact that 7 ^||c°||a,ft, + 2($(c°), 1) < 1. To prove this bound, observe 
that 

7 ^||c°|U,, + 2(<i>(c°), 1) < 7 ^||c°||?,, + 1 + ||c°||i4(o) + kir 

where we have used the definition (2) of the free-energy <i> in the first line followed by 
the discrete Friedrichs’ inequality with r = 4, 2 in the second line and the first bound 
on the initial datum in (46) below to conclude. 
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(ii) We next prove that 


(35) 


N 

E 




The discrete Agmon’s inequality (28) followed by the first inequality in (22) yields 


N 


N 




n=l 


n=l 


N 


n=l 


The first factor is < 1 owing to (30). Thus, to prove (35), it suffices to show that 
also the second factor is < 1. Using the definition (27) of followed by (24b) with 
we infer that 


(36) 


I^WL'^hcIWIh = = ($'(c)(), L^c)() - «,L^c)(). 


T k ri 


Using again (27) for the second term in the right-hand side of (36) followed by the 
Cauchy-Schwarz and Young’s inequalities, we obtain 




Hence, since \wl\o,h < ^ W^WaM, 


The fact that 2^1 ~ ^ then follows multiplying the above inequality by 

r, summing over 1 ^ n ^ A^, using (30) to bound the second and third term in the 
right-hand side, and observing that 


(37) ||‘i>'(cDP<l|c: 


h\\L<i(Q.) + 2 || 


,n||4 

'h\\L'^(Q.) 


+ c: 


.rail 2 


< 


c”ll® 


+ C 


,77 II 4 


^h\\l,h 


+ c: 




2 

1,/t 


< 1 , 


where we have used the definition (2) to obtain the first bound, Friedrichs’ inequal¬ 
ity (18) with r = 6,4, 2 to obtain the second bound, and (30) together with the first 
inequality in (22) to conclude. 

(hi) We proceed by proving that 


(38) max Kp4-y2|’^||^^,ra||2<i^ 
l^n^N 

n=l 

Let := — q^Acg). Recalling (25), satisfies 

(39) {wl,'il}h) = {^'{cl),'tph) + l'^ah{dL,±f) 

For any 1 ^ n ^ A^, subtracting from (24b) at time step n (24b) at time step (n — 1) 
if n > 1 or (39) if n = 1, and selecting = w)) as a test function in the resulting 
equation, it is inferred that 

« - = rjW{Stcl,wl) + - $'(cr!),<). 
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Using (24a) with to infer = -Tj'^\\Stcl\\'^, we get 

(40) « - = ($'(c)() - $'(cri),0. 

From the fact that 


(41) 


<i)'(r) — 4)'(s) = (r^ + rs + — l)(r — s), 


followed by the Cauchy-Schwarz and Young’s inequalities, we infer 

T-'V^ tU" 

(42) i(4>'(c)() - $'(cri),oi ^ -^\\s,cir + —iKir, 

with C" := C{1 + ||c))||^a)(f 2 ) + l|c/(~^|li<x.(o)) for a real number C > 0 independent of 
h and r. Using (33) for the first term in the left-hand side of (40) together with (42) 
for the right-hand side, we get 

(43) IKP -UK - K-'P + ^ rC'^IIKP + 

Summing (43) over 1 ^ n ^ IV, observing that, thanks to (35) and the second bound 
in (46) below, we can have tC” < 1 for all 1 ^ n ^ provided that we choose r small 
enough, and using the discrete Gronwall’s inequality (29) (with <5 = t, a" = 
jfi ^ ^ (jn Q ^ ||^ 0 || 2 ^^ estimate (38) follows if we can bound 

To this end, recalling the definition of and using the Cauchy-Schwarz 
inequality, one has 

lk°P = -7 ^(Aco,K s: (||4’'(c^)|| + i'^WcoWh^q)) lk°||, 

and the conclusion follows from the regularity of cq noting the first bound in (46) 
below and estimating the first term in parentheses as in (37). 

(iv) We conclude by proving the bound 

(44) max (\K\\L’^(n) + \\Ltcl\\l,h) ^ 1- 

l^n^N \ / 


Using the Cauchy-Schwarz and Young’s inequalities to bound the right-hand side 
of (36) followed by (18) with r = 6,4,2 and the first inequality in (22), we obtain, for 
all 1 ^ n ^ A^, 


(45) 


7^l|i^c^llo.K7-Ml|4>'(c)()P+||Kp) 


^ {KWUn) + KWUn) + \Kf) + Kir 

~ (l|c?(||a,?i + l|c))||a,^ + -I- IIKir ^ 


where we have concluded using (30) multiple times for the terms in parentheses 
and (38) for ||w^p. Using the discrete Agmon’s inequality (28) followed by Young’s 
inequality and the first inequality in (22), we infer 


max l|c”||ia>(n) < max (||c)(|U,/^ -f ||tKIIo./«) 

l^n^N l^n^N \ / 


< 1 , 


where the conclusion follows using (30) for the first addend in the argument of the 
maximum and (45) for the second. □ 



12 


F. CHAVE, D. A. DI PIETRO, F. MARCHE, F. PIGEONNEAU 


Proposition 9 (Bounds for c°). Let c° e be defined by (25) from an initial 
datum Co e n L‘^{Vl) such that = 0 on dVl. It holds, with real number 

C > 0 independent of h, 

(46) l|c/ll|i,/i + ||c°||_Laj(o) ^ C'||co||_f/ 2 (Q). 

Proof. To prove the first bound in (46), let in (25) and use the first inequality 

in (22), the Cauchy-Schwarz inequality and the discrete Poincare’s inequality (18) 
with r = 2 to infer 


llc^ll?,. ^ a,(c°,c°) = -(Aco,c°) < ||Aco||||c°|| < ||cob.(a)||c°||i.,. 

To prove the second bound in (46), we start by noticing that, using the definition (27) 
of LX with = -IXqI, 

llA^c°||o_^ = —ahiXh, fXXh) = {/Lco,L^fil) ^ ||co||/i-2(Q)||I/^c)(||, 

hence ||T^C/(||o,/i < ||co||r 2 (o). Combining the discrete Agmon’s inequality (28) with 
the latter inequality and the first bound in (46), one gets 

||c°||lco(o) ^ ||c°||i_^||L^C°||^;^ < ||co||_H-2(q), 

and the desired result follows. □ 


5. Error analysis. In this section we carry out the error analysis of the method (24). 

5.1. Error equations. Our goal is to estimate the difference between the dis¬ 
crete solution obtained solving (24) and the projections of the exact solution such 
that, for all 1 ^ n ^ N, = fXw^, while, for all 0 ^ n ^ A^, c)( e UXfi solves 

a/i(cft,^^) =-(Ac”,(/3,i) yPf^eUX, 

and (cj), 1) = 0. We define, for all 1 ^ n ^ A^, the errors 

(47) elh--=cl-cl, el,H--=wl-wl 

By definition (25), c° = c°, which prompts us to set e° ^ := 0. Using Poincare’s 
inequality (18) with r = 2 and the consistency (23) of Oh, the following estimate is 
readily inferred: For all 0 ^ n ^ A^, assuming the additional regularity c" G iJ*’+^(n), 


(48) < K-I^c-lli,, <h'=+i||c"||^...(0). 

Remark 10 (Improved L^-estimate). We notice, in passing, that, using elliptic regu¬ 
larity (which holds since 11 is convex, cf, e.g., [25]), one can improve this result and 
show that ||c^ — < /i^+^||c"||jjfc+ 2 (f 2 ). 

Recalling (24), for all 1 ^ n ^ A^, the error (e"^,e” ^) G g ^ UX solves 

(49a) iph) + = ^(^^) ^ 

(49b) {el^u.^h) = {^'{cD - l>'{c^)M + ^ 
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where, in (49a), we have defined the consistency error 
(50) £{'£^) ■■=-{5t^,Lph) - 

while in (49b) we have combined the definitions of wJJ and with (lb) to infer 

{wl,'il)h) - = (m” + Ac”, ■!/'/«) = ($'(c”),'i/'/i)- 


5.2. Error estimate. 

Theorem 11 (Error estimate). Suppose that the assumptions of Lemma 8 hold true. 
Let {c,w) denote the solution to (1), for which we assume the following additional 
regularity: 

(51) c e C2([0, tF]-,L\n)) n ^^([O, tp]; ff"+"(f^)), w e C°([0, If]; . 


Then, the following estimate holds for the errors defined by (47); 


(52) 


N 


max 

l^n^N 


^c,h\\ 


+ 


1] T\K,h\\ 


71=1 


^ C{h^+^ + r). 


with real number C > 0 independent of h and r. 

Proof. Let 1 ^ n < TV. Subtracting (49b) with from (49a) with 

we obtain 


(53) lie”,,111, + 7"a,(el„ (5tel,) = £:(e”,,) + ($'(c”) - $'(c)(), 5*6^,) := Ti + 


We proceed to bound the terms in the right-hand side. 


(i) Bound for Ti. Let e U'f. Adding to (50) the quantity 

(dtc” - Au;”, ph) + (^iTT^+^c” - Stc^, pu) = 0, 

(use (la) to prove that the first addend is 0 and the definition of the L^-orthogonal 
projector to prove that the second is also 0), we can decompose as follows: 

£{p^) = (dtc” - + (5t(7r^+ic” - cff),ph) - [ah{wl,p^) + (Aw^pn)) 

■= ‘^1,1 + ^1.2 + 2^1,3- 
For the first term, we have 

(54) ITpil sS ||dtc”-(5tc”||||(p,|| <'r||c||c 2 ([o,tp];L 2 (n))ll^;,l|i.h ^ rllpj^t, 

where we have used the Cauchy-Schwarz inequality, a classical estimate based on 
Taylor’s remainder, Poincare’s inequality (18) with r = 2, and we have concluded 
using the regularity (51) for c. For the second term, on the other hand, using the 
Cauchy-Schwarz inequality followed by (48) together with the C'^-stability of the 
backward differencing operator (3), Poincare’s inequality, and the regularity (51) for 
c, we readily obtain 

(55) < m7T^^+^c--c-J\\\\p^\\ < /r"+l|c”||ci([o.*p];H'=+2(o))||^,|| < 
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Finally, recalling the consistency properties (23) of Uh, we get for the last term 


(56) 






Collecting the bounds (54)-(56), it is inferred that 

(57) $ := sup ^ 

so that, for any real e > 0, denoting by Cg > 0 a real depending on e but not on h or 
r, and using the second inequality in (22) to bound ||ej); ^ ||e(), ^||a,/i, 

(58) |Ti| < $||e-,,||i,. < (h'^+i + r)||e(L,J|i., ^ + r)^ + 


(ii) Bound for T 2 . Set, for the sake of brevity, Q" := ^'(c))) — $'(c"), and define the 
DOF vector e U_'f such that 

(59) 

zy = 4+ig" yreTh, ^f = 4{Q”}f ZF = ^'fzT^ VF’gJ-^ 

where {-If denotes the usual average operator such that, for any function ip admitting 
a possibly two-valued trace on F e Fti fi J-T 2 , W}f '■= ^{‘P\Ti + while, for a 

boundary face F e Tp denotes the unique element in Th such that F e Fpp- 
We have, using the definition of followed by (49a) with = z^^, (57), and the 
second inequality in (22), 

(60) T 2 = {zh,Stelf,) = £{Zf,)-ahiel^i„Zf,) < ($ + \\el^i,\\a,h) khl|i,/i- 
By Proposition 14 below, 

(61) \\zJi,H<\\elJa,h + h>^+\ 

hence, for any real e > 0, denoting by Ce > 0 a real number depending on e but not 
on h or r, and recalling the bound (57) for $, 

(62) |‘I 2 | < (l|e"/illa,h + + r)^) + e||e” 


(iii) Conclusion. Using (58) and (62) with e = | to bound the right-hand side of (53), 
it is inferred 


hl,h\\l,h + ^ (h'^+Vr)^ 


2 

a,h‘ 


Multiplying by r, summing over 1 ^ n ^ A^, using (33) for the second term in the 
left-hand side, and recalling that, by definition, e(!^ = 0, we get 


N 

'y^hc.h\\l,h+ 2 ^l|e”,/. 

n=l 


N 


< 


T.Cr\\e:jlH + C{h'^^^+rr, 


with C > 0 independent of h and r. The error estimate (52) then follows from an 
application of the discrete Gronwall’s inequality (29) with (5 = r, a" = 7 ^||h 


6 " = 


p X” = C”: and G = C'(/i^+^ 4- r)^ assuming r small enough. 


c,h \\a,h ’ 

□ 
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Remark 12 (BDF2 time discretization). In Section 6, we have also used a BDF2 
scheme to march in time, which corresponds to the backward differencing operator 

4 V:- - 2 ;-. 

used in place of (3). The analysis is essentially analogous to the backward Euler 
scheme, the main difference being that formula (33) is replaced by 

2x{3x — 4y + z) = x^ — + {2x — y)^ — {2y — z)^ + {x — 2y + z)^. 

As a result, the right-hand side of (52) scales as [hf'^^ + r^) instead of + r). 

To prove the bound (61), we need discrete counterparts of the following Gagliardo- 
Nirenberg-Poincare’s inequalities valid for p e [2, + 00 ) if d = 2, p e [2,6] if d = 3, and 
alive H^{n) nLlin): 

( 63 ) ^ ^ ^ ^ ’ 

where the first bound follows from [1, Theorem 3] and the second from Poincare’s 
inequality. The proof of the following Lemma will be given in Appendix A. 

Lemma 13 (Discrete Gagliardo-Nirenberg-Poincare’s inequalities). Under the as¬ 
sumptions of Lemma 6, it holds for p e [2, +cx)) if d = 2, p e [2,6] if d = 3 with C > 0 
independent of h and a defined as in (63), 

(64) e Ulo, W^HVhhHny ^ 

Proposition 14 (Bound on ||z^||i_/i). With Zf^ defined as in (59), the bound (61) 
holds. 

Proof. Recalling the definition (15) of the ll-Hi /,-norm, one has 

(65) 

TeTh FeJ^Tr^J^f 


(i) Bound for Ti. Using the id^-stability (7) of , formula (41) to infer Q” = 
— c") with := (c)))^ -f -I- (c")^ — 1, the triangle and Holder inequalities, 
we get, for all T e Th, 

12:11 < ||V/,g"|l < - c")|| + 11(c)) - c")V/,g"|| 

^ (l|c))|lico(n) + l|c"||ioo(o) + 1) \\Vh{cI - c'‘)|| 

+ l|c/) — c"'||L8(n) (||c))||ioo(Q) -I- ||c"’||ioo(Q)) (||V/ic))||/,3(n)d -f ||Vc"||i3(o)d) . 

Noting the a priori bound (44) and the regularity assumption (51), both ||c))||icx)(Q) and 
||c"||ioo(n) are < 1. Additionally, by the continuous Gagliardo-Nirenberg-Poincare’s 
inequality (63) with p = 3 and the regularity assumption (51), one has with a = 
1 / 2 -I-‘yi 2 , II Vc" 11^3(0)d < |c"|^“q)||c”||^ 2 (q) ^ 1- Similarly, the discrete Gagliardo- 
Nirenberg-Poincare’s inequality (64) with p = 3 combined with the a priori bounds (30) 
and (44) yields || V/,c))||i3(n)d < ||c))||);;“||L^c))||));, < 1. Then, inserting ±(8))-7r^+ic”) 
and using the triangle inequality, 

|Ti| < {WuelA + VlnWa)) + (l|V/,(8)) ~ + \\cf, - ^r'c"|Ua(o)) 

+ (II V/,( 7 r^+ic" - c")II + IItt^+Ic" - c”||ia(o)) := Ti,! + 2:1,2 + 2:1.3. 
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Using the discrete Friedrichs’ inequality (18) with r = 6 together with the defini¬ 
tion (15) of the ||•||l.;i-norm and the first inequality in (22), it is readily inferred 
that Ti,i < Again the Friedrichs’ inequality (18) with r = 6 followed 

by the approximation properties (48) of and the regularity (51) yields T 2,2 ^ 
||c”||j:/)=+ 2 (Q) < Finally, using the approximation properties (8) of we 

have Ti ^3 < h*’+^(||c"||//fc+ 2 (Q) -I- ||c"||^yfe+i,6(o)) < where we have used the fact 

that i/^+^(U) c VF^+^’®(U) for all /c > 0 and d e {2,3} on domains satisfying the 
cone condition (cf. [2, Theorem 4.12]). Gathering the previous bounds, we conclude 
that 

(67) \T,\<\\el4a,H + h>^+\ 

(ii) Bound for T 2 - For all interface F e Ft^ n we denote by [-{f the usual 
jump operator such that, for every function ip with a possibly two-valued trace on 
F, \_p\f '■= ^\Ti ~ V\T 2 (tFe orientation is irrelevant). Let an element T e Th and an 
interface face F e Ft n Ft+ be fixed. Using the L^-stability of tt^, inserting ±Qt 
( with Qt Q”|t)i and using the triangle inequality it holds, 

II 4 ({Q”}f - < II{Q”}f - f 44 ?IIf 

(68) ^ll[Q"]F||F + ||Q^-7r44TllF 

< ||[Q"]f||f + 4||VQJ||t, 

where we have used (8) for the last term. Let us bound the first term in the right-hand 
side. Observing that [<i)'(c")]F = 0 and recalling (41), it is inferred 

|[Q"]f| = |[$'(c}()]f| < IKJfI (|ct|^ + |cT||cr.| + |CT.|^ + 1). 

Using this relation, and noticing the a priori bound (44), we get 

I|[Q"]f||f ^ (l|c}}||icc.(o) +1^ IIKIfIIf ^ IIKIfIIf = IK - c”]f||f, 

where the conclusion follows observing that c" has zero jumps across interfaces. In¬ 
serting + [c}} — f4^c”]f inside the norm and using the triangle inequality, we obtain 

(69) ||[Q"]f||f < IK - F^fWf + m - f4'c”]f||f + - c"]f||f. 

Define on H^{Th) the jump seminorm jujj := Z;F€Fj^ ^f^II[^]f|If' Let us prove that 

(70) e Ut \vh\j < \\vji,h ^ \\vh\\a,h- 

Inserting ±(Fp[uh]F — vf) and using the triangle inequality, it is inferred that 

^ Zl Z ^FH\\^T-TrFVT\\F+hUvT-VF)\\F) <\\VhVhV + \v^\lh, 

F€F‘ Fen- 


where we have used (9) followed by the discrete trace inequality (4) and the fact that 
card(J^F) < 1 by mesh regularity for the first term, and the definition (15) of the 
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Fig. 1: Mesh families for the numerical tests 


Mi,?i-seminorm for the second term. This proves the first bound in (70). The second 
bound follows from ( 22 ). 

_ 1 

Multiplying ( 68 ) by , squaring, summing over F e then over T B Th, using 

mesh regularity to infer that card(J^T) is bounded uniformly in h, and noticing (69) 
yields 


(71) 


^ l|V,g"p + K - + |c)( - - c"|2 

^ II + \\cl- Ltcllln + - c"|2 

< l|V,g"p + ||e^,J|2 , + (h^-+i||c"||^...(n))' , 


where we have used (70) to pass to the second line and the approximation proper¬ 
ties (48) of c)( and ( 8 ) of to conclude. Proceeding as in point (i) to bound the 
first term in the right-hand side of (71), and recalling the regularity assumptions (51) 
on c, we conclude 


(72) 


1^2 I < lle^.JIa,^ 


+ h'^+\ 


(iii) Conclusion. Using (67) and (72) in (65), the estimate (61) follows. □ 

Remark 15 (Polynomial degree for element DOFs). The use of polynomials of degree 
{k + 1) (instead of k) as elements DOFs in the discrete space (13) is required to infer 
an estimate of order in ( 66 ) and for the last term in (71). 

6. Numerical results. In this section we provide numerical evidence to confirm 
the theoretical results. 


6.1. Convergence. We start by a non-physical numerical test that demon¬ 
strates the orders of convergence achieved by our method. We solve the Cahn-Hilliard 
problem (49) on the unit square D = (0,1)^ with tp = 1, order-parameter 

(73) c{x,t) = tcos{TTXi) cos{Trx 2 ), 

and chemical potential w inferred from c according to (lb). The right-hand side of (la) 
is also modified by introducing a nonzero source in accordance with the expression of 
c. The interface parameter 7 is taken equal to 1. 

We consider the triangular, Cartesian, and (predominantly) hexagonal mesh families 
of Figure 1. The two former mesh families were introduced in the FVCA5 bench¬ 
mark [27], whereas the latter was introduced in [20]. To march in time, we use the 
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Fig. 2: Energy-errors at final time vs. h. From left to right: triangular, Cartesian 
and (predominantly) hexagonal mesh families; cf. Figure 1. 


implicit Euler scheme. Since the order-parameter is linear in time, only the spatial 
component of the discretization error is nonzero and the choice of the time step is 
irrelevant. The energy errors and at final time are 

depicted in Figure 2. For all mesh families, the convergence rate is (A: -(- 1), in accor¬ 
dance with Theorem 11. For the sake of completeness, we also display in Figure 3 the 
L^-errors ||c)) — and ||w)( — , for which an optimal convergence rate 

of (fc -f 2) is observed. 


6.2. Evolution of an elliptic and a cross-shaped interfaces. The numerical 
examples of this section consist in tracking the evolution of initial data corresponding, 
respectively, to an elliptic and a cross-shaped interface between phases. For the elliptic 
interface test case of Figure 4, the initial datum is 

fo.95 if81(a;i-0.5)^ + 9(a;2-0.5)^ < 1, 

cnia;) = { 

I —0.95 otherwhise. 


For the cross-shaped interface test case of Figure 5, we take 

if 5 (|(a::2-0.5) - |(a:i-0.5)| + ||(xi-0.5) + (x2-0.5)|) < 1 


co{x) = < 


0.95 


-0.95 


or 5(|(xi—0.5) 
otherwhise. 


g V J- j I ' I 5 

|(x2-0.5)| + '^ 


(x2-0.5) + (xi-0.5)|) < 1, 


In both cases, the space domain is the unit square = (0,1)^, and the interface 
parameter 7 is taken to be 1 • 10“^. We use a 64 x 64 uniform Cartesian mesh and 
k = \ with time step r = 7^/10. 
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Fig. 3: L^-errors at final time vs. h. From left to right: triangular, Cartesian and 
(predominantly) hexagonal mesh families; cf. Figure 1. 


In the test case of Figure 4, we observe evolution of the elliptic interface towards a 
circular interface and, as expected, mass is well preserved (+0.5% with respect to the 
initial ellipse). Similar considerations hold for the cross-shaped test case of Figure 5, 
which has the additional difficulty of presenting sharp corners. 


6.3. Spinodal decomposition. Spinodal decomposition can be observed when 
a binary alloy is heated to a high temperature for a certain time and then abruptly 
cooled. As a result, phases are separated in well-defined high concentration areas. In 
Figure 6 , we display the numerical solutions obtained on a 128 x 128 uniform Cartesian 
mesh for fc = 0 and on a uniform 64 x 64 Cartesian mesh for k = 1. In both cases, we 
use the same initial conditions taking random values between -1 and 1 on a 32 x 32 
uniform Cartesian partition of the domain. The interface parameter is 7 = 1/100, 
and we take r = 7 ^/ 10 . For fc = 0, the time discretisation is based on the Backward 
Euler scheme while, for fc = 1, we use the BDF2 formula to make sure that the spatial 
and temporal error contributions are equilibrated; cf. Remark 12. 

The separation of the two components into two distinct phases happens over a very 
small time; see two leftmost panels of Figure 6 corresponding to times 0 and 5 • 10“®, 
respectively. Later, the phases gather increasingly slowly until the interfaces develop 
a constant curvature; see the two rightmost panels of Figure 6 , corresponding to times 
1.25 • 10“^ and 3.6 • 10“^, respectively. At the latest stages, we can observe that the 
solution exhibits a (small) dependence on the mesh and/or the polynomial degree, 
and the high-concentration regions in Figures 6 a and 6 b are highly superposable but 
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Fig. 4: Evolution of an elliptic interface (left to right, top to bottom). Displayed 
times are 0,3- 10“^ , 0.3, 1. 




Fig. 5: Evolution of a cross-shaped interface (left to right, top to bottom). Displayed 
times are 0, 5 • lO-^, 1 • lO’^, 8.17 • IQ-^. 


not identical. 

Appendix A. Proofs of discrete functional analysis results. 

This section contains the proofs of Lemmas 6 and 13 preceeded by the required pre¬ 
liminary technical results. 
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Proposition 16 (Estimates for L^). 

(74) G C/t 

(75) Vi;, g 


Assuming mesh quasi-uniformity (10), it holds 

\\lXvJo,h ^ 

WLkUhWH-^Q) ^ 


Proof, (i) Proof of (74). Let e UX- Making z, = —lXv.h the definition (27) of 
fX, we have 


W^X.ldhWo.h 


= -ah{: 


’-hJ If-hlth) 


< 


^hWiALhUhhh 


< 


ithWiM 


IIL 


h'l’-hWo.h, 


where we have used the continuity of Oh expressed by the second inequality in (22) 
followed by the fact that, for all z, e UX, |U;i||i,/i ^ h~^\\zXo,h- This inequality 
follows from the definition (15) of the H-Hi^ii-norm using the inverse inequality (5) to 
bound the first term and recalling mesh quasi-uniformity (10). 

(ii) Proof of (75). Let u, e UXq. Observing that has zero-average on fl (cf. 

Remark 5), we have 


(76) \\LlvXH~^n)= sup (L^u,,^?). 


Let now := Xp- Using the fact that e P^'+^(7)i) followed by the defini¬ 

tions (27) of fX and (26) of (-jOo.Iij one has 

{Llvi,,ip) = {L'fvh,TrXX) = -so.ft(T^L/i,^,) - a/i(llft,^,)- 
Hence, using the Cauchy-Schwarz inequality we get 

\{LhV^,ip)\ < \L^vXo,h\pX°A + IIV/iII 1.1*11^,IIi./i 

^ h~^hh\\i,hh\ipXi,h + I|t/^I|i./«II^,I|i./i 

^ llll/illi.ft||^,lli./i ^ llll/.lli.ft|IV5||Ri(n), 

where we have used the second inequality in (22) in the first line, (74) together with 
the fact that |z,|o.ii < for all z, e iX to pass to the second line, and the 

iV^-stability (17) of X fo conclude. To obtain (75), plug the above estimate into the 
right-hand side of (76). □ 

We introduce the continuous Green’s function Q : Lo(U) H^{n) n Lo(U) such that, 
for all ip e Lq(H), 

(Vt/(/ 3 , Vu) = ((/;, u) VuGi 7 ^(r 2 ). 

Owing to elliptic regularity (which holds since H is convex), we have Qip e Its 

discrete counterpart X • U^,o ^ U defined such that, for all e UX,o^ 

(77) ^h{G_f^‘P^T ^h) ~ {Pf^i^h)o,h Vz, G G, Q, 

with inner product {■,-)o^h defined by (26). We will denote by G^Vf^ (no underline) 
the broken polynomial function in P^+^(7)i) obtained from element DOFs in Gh-h- 
We next show that —is the inverse of fX restricted to UX,o UXg- Let u, G UXg- 
Using (77) with ip^ = fXv-h followed by (27), it is inferred, for all z, G UX,o^ 

ahiX^^hV-h^^h) = iLtvh^Zh)o.h =-ah{v^,zX aXllh + OX^^hV-h^^h) = 0- 
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Therefore, since + G'^LhV.h) ^ U-h.o ah is coercive in Uh^Q (cf. (22) and Propo¬ 
sition 2), we conclude 

( 78 ) Vh + Q.h^hV.h=Q- 

Proposition 17 (Estimates for Q ^). It holds, for all Vh e U_h,o: 

(79) \\Q!hlIh - LhG'^h\\l,h ^ h {\Vh\oM + \\GVh\\H2{Q)) ^ h\\]lh\\o,h- 
Moreover, using elliptic regularity, we have 

(80) WOhUh - TTh'^^GvhW < h^ {\vh\o,h + WGvhWH^n)) ^ h'^hh\\o,h- 
Proof. Let e Uhfi¬ 
ll) Proof of (79). For all Zf, e U'f q we have, using the definition (77) of Gh-h 
subtracting the quantity {vt + AGvh,Zh) = 0, 

(81) ahiGl-Hh - IhGvh,Zffl = {vi^,Zff)o,h - {vh, Zh) -ahUlGvh, Zffj - {AGvh, Zh) ■ 

' - , -''-.-' 

'Xi X2 

Recalling the definition (26) of the inner product (•, one has 

(82) |Ti| = \so,h{vh,Zh)\ < \vh\e),h\zh\o,h < h\vh\o,h\zh\i,h- 

On the other hand, the consistency property (23) of the bilinear form ah readily yields 


(83) IT 2 I < h\\Gvh\\H^(n)\\zi^\\i,h- 

Making = G\V-h~^X^'*^h in (81), and using the coercivity of ah expressed by the first 
inequality in (22) followed by the bounds (82)-(83), the first bound in (79) follows. To 
prove the second bound in (79), use elliptic regularity to estimate \Gvh\H'^(n) ~ l|i^h|| 
and recall the definition of the ||•||o,Ii-norm. 

(ii) Proof of (80). We follow the ideas of [19, Theorem 10] and [18, Theorem 11], to 
which we refer for further details. Set, for the sake of brevity, := Gh^h ~ 
and let z := Gph- By elliptic regularity, z e H^{Q) and 1 |z1|//2(q) < l|y)?il|. Observing 
that —Az = iph, letting z^ := /^z, and using the definition (77) of we have 

(84) WphW^ = -{Az,iph) - ah{(p,,Zff) + {vh,Zh) - ah{,lXGvh,z,^) + sofi{v^,z_^). 

^ ^ ^ '-.-' '-.-' 

Ti “^2 T 3 

Using the consistency (23) of ah, it is readily inferred for the first term 

(85) |Ti| ^ /il|zl|//2(f2)||^^||i_,j < h'^ {\vh\o,h + \\Gvh\\H2(n)) ||<7’ii||, 

where we have used elliptic regularity to infer ||z:||// 2 (q) < ||(p;i|| and (79) to bound 
For the second term, upon observing that (u/i, z/j) = —{AGvh,z) = {'VGvh,'^z) 
since, by definition of, —AGvh = Vh ^ P^+^(7)i) and Zh = Tr’fG'^z, recalling the def¬ 
inition (21) of the bilinear form ah and using the orthogonality property (20) of 
(Pt^^ ° Lt)^ '"^6 have 

^2 = (V(p^+lj^0U/, - GVh),V{p!f+^Zh - z))t + Sfih{dGVh,Zh). 

Ten 
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By the approximation properties of ° Lt) of ^ bounding ||z||// 2 (q) 

and ||<p^||i,/t as before, we have 

(86) Y^2\ ^ {\'^h\o,h + \\Guh\\H^(n)) IPhl- 

Finally, for the last term, we write 

(87) l^al < \vh\o,h\zh\o,h ^ \vh\o,hh?\\z\\H^n) ^ h^\vh\o,hWh\\, 

where we have used the Cauchy-Schwarz inequality in the first bound, the approxima¬ 
tion properties (8) of in the second bound, and elliptic regularity to conclude. 
Using (85)-(87) to estimate the right-hand side of (84) the first inequality in (80) 
follows. Using elliptic regularity to further bound \Gvh\H'^{p,) ^ ||uft|| and recalling 
the definition of the ||-lo./t-norm yields the second inequality in (80). □ 

Remark 18 (Choice of So,/i). The choice (26) for the stabilisation bilinear form 
so,h is crucial to have the right-hand side of (87) scaling as h^. Penalizing the full 
difference {vp — vp) instead of the lowest-order part TTp{vF — vp) would have lead to 
a right-hand side only scaling as h. 

We are now ready to prove Lemma 6. 

Proof of Lemma 6. Let e U^ q and set := Recalling that, owing to (78), 

Vh = it is inferred using the triangle inequality, 

(88) ||^^?i||L'»(n) < ~ '■= + ‘^ 2 - 

The L”-stability of (cf. (7)) followed by the continuous Agmon’s inequality 

readily yields for the first term 

(89) Ti < \\G‘Ph\\L’^(n) ^ ll07’?ill|fi(n)ll0'^ft||^2(Q)- 

Using a standard regularity shift (cf., e.g., [25]), recalling that ph = and using 

the i7“^-bound (75) for L^Vf^, we have 

(90) ^ ll7’/i||R-i(n) ^ WGphWn^in) ^ \\Ph\\ = \\LhV}^\\, 

which plugged into (89) yields 

(91) 


For the second term we have, on the other hand. 



where we have used the global inverse inequality (12) with p = 2 to obtain the first 
bound, the estimate (80) to obtain the second, (74) to obtain the third, and the fact 
that d ^ 3 together with h ^ Lq < 1 (with /in diameter of LI) to conclude. The 
conclusion follows plugging (91) and (92) into (88). □ 
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Remark 19 (Discrete Agmon’s inequality in dimension d = 2). When d = 2, we have 
the following sharper form for the discrete Agmon’s inequality: 


(93) 




Hh 





To obtain (93), the following modifications are required in the above proof: (i) The 
term^i is bounded as Ti < < ||'a/i|| ® || 2 ^ where we have used 

Vfi = —Gq^h (cf. (78)^ for the first factor and (90) for the second; (ii) The third line 

of (92) becomes ‘I 2 < {h\\vk\\i,h)^\\LkVk\\lk ^ hh\\o.h\\^^hTh\\lh> we have used 

the inverse inequality (5) and mesh quasi-uniformity to bound the first factor. 


We next prove the discrete Gagliardo-Nirenberg-Poincare’s inequality of Lemma 13. 
Proof of Lemma 13. Using the same notation as in the proof of Lemma 6, we have 


\\^hVh\\Lp{a)‘‘ ^ hT^h^^G'Ph\\Lp{n)‘‘ + h{Gh'f_k~ '^h^^Gq:h)\\Lp{ny Ti + T2- 

For the first term, we use the VF^’^-stability of (cf. (7)) followed by the continuous 
Gagliardo-Nirenberg-Poincare’s inequality (63), and (90) to infer 

For the second term, on the other hand, we have 




< 


h<-p--^)- tMWfTKPk - LhGThWl 
<h-+<^-^^){h\\Ltvk\\o,Hy 

< /l“+<|-5)|( 


WLLhUhWo.h 




JjI j 


iil-ai 


where we have used the global reverse Lebesgue inequality (11) in the first line, the 
definition (15) of the ||•||l.ft,-norm to pass to the second line, the estimate (79) to pass 
to the third line, and (74) to pass to the fourth line. To obtain the second inequality in 
the fourth line, we observe that, recalling the definition (63) of a and the assumptions 
on p, it holds for the exponent of /i. 


.7(1 1 

a + d\ -- 




and, since ^ /iq < 1, the conclusion follows. □ 

Remark 20 (Validity of the discrete Agmon’s and Gagliardo-Niremberg-Poincare’s 
inequalities). At the discrete level, the fact that the discrete Agmon’s inequality (28) is 
valid only up to d = 3 and that the Gagliardo-Nirenberg-Poincare’s inequalities (64) 
are valid only for p e [2, +cx)) if d = 2, p e [2,6] if d = 3 is reflected by the need 
to have nonnegative powers of h in the estimates of the terms %2 to conclude in the 
corresponding proofs. 
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(a) 128 X 128 uniform Cartesian mesh, A: = 0, BE 




(b) 64 X 64 uniform Cartesian mesh, fc = 1, BDF2 


Fig. 6: Spinoidal decomposition (left to right, top to bottom). In both cases, the 
same random initial condition is used. Displayed times are 0, 5 • 10“®, 1.25 • 10“^, 
3.6 • 10-2. 










